Can noncommutative effects account for the present speed up of the cosmic 

expansion? 
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In this paper we investigate to which extent noncommutativity, a intrinsically quantum property, 
may influence the Friedmann-Robertson- Walker cosmological dynamics at late times/large scales. 
To our purpose it will be enough to explore the asymptotic properties of the cosmological model in 
the phase space. Our recipe to build noncommutativity into our model is based in the approach of 
reference [Phys. Rev. Lett. 88 (2002) 161301], and can be summarized in the following steps: i) the 
Hamiltonian is derived from the Einstein-Hilbert action (plus a self-interacting scalar field action) for 
a Friedmann-Robertson- Walker spacetime with flat spatial sections, ii) canonical quantization recipe 
is applied, i. e., the minisuperspace variables are promoted to operators, and the WDW equation is 
written in terms of these variables, iii) noncommutativity in the minisuperspace is achieved through 
the replacement of the standard product of functions by the Moyal star product in the WDW 
equation, and, finally, iv) semi-classical cosmological equations are obtained by means of the WKB 
approximation applied to the (equivalent) modified Hamilton-Jacobi equation. We demonstrate, 
indeed, that noncommutative effects of the kind considered here, can be those responsible for the 
present speed up of the cosmic expansion. 
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I. INTRODUCTION 



Noncommutativity of spacetime coordinates - an old 
idea dated back to 1947 [l[ - is the simplest expected 
modification to quantum field theory. Noncommutativ- 
ity - the central mathematical concept in quantum me- 
chanics - expresses uncertainty in the simultaneous mea- 
surement of any pair of conjugate variables, such as po- 
sition and momentum. In the presence of a strong mag- 
netic field noncommutativity arises 0], even in a clas- 
sical context. More recently, noncommutativity has re- 
ceived increased interest in connection with developments 
in string theory. Attempts to connect M(atrix)-string 
theory to cosmology on the brane Q have shown that 
noncommutativity arises in the former theory. 

There are several approaches in the literature to build 
noncommutativity into field theories. One of the au- 
thors and collaborators have explored several of these 
approaches 0413 • ^ n some of these formalisms the as- 
sumption of noncommutativity among the spacetime co- 
ordinates has a consequence that the fields present do 
not commute themselves. This is the particular case of 
the Seiberg-Witten map I13ll and its generalization by 
Wess and collaborators [l4j, in which noncommutative 
fields are obtained as an (infinite) expansion of the usual 
commutative fields in the noncommutative parameter. 



Space-time noncommutativity has as a consequence a 
different product of functions (the Moyal/star product), 
which induces noncommutativity among the fields. 

In cosmological settings an already well explored way 
to include the effects of noncommutativity is given in 
reference [H , where deformation of spacetime itself is re- 
placed by noncommutativity in minisuperspace instead 
(see also references [HI, [l6[ where similar approaches are 
applied). 1 The latter proposal is inspired in various re- 
lated results in the literature, among which we can cite 
the above mentioned case of the Seiberg-Witten map 
[1> EH, EH ■ In the case of quantum cosmology, the min- 
isuperspace variables play the role of the coordinates in 
configuration space. 2 Thus, as stated in it seems rea- 
sonable to propose a kind of noncommutativity among 
these specific gravitational variables. The noncommu- 
tative proposal there is formulated in terms of models 
with a finite number of degrees of freedom, where the 
Wheeler-DeWitt (WDW) equation resembles a Klein- 
Gordon equation, this time in terms of the minisuper- 
space variables. One can then apply the same procedure 
as in noncommutative quantum mechanics to end u p w ith 
a noncommutative version of quantum cosmology [2l| . 

In the present paper we shall investigate the possible 
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1 For further consideration of the approach of Q in the context 
of extra-dimensional cosmology see, for instance, Ref. |17H and 
references therein. 

2 See, for instance, references [l9, 20], for an alternative path in- 
tegral approach to quantum cosmology. 
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effect of noncommutativity in the late-time/large scale 
cosmological dynamics. In reference [22| the effects of 
noncommutativity were investigated within a dilatonic 
cosmological model for an exponential dilaton potential. 
The existence of such noncommutativity results in a de- 
formed Poisson algebra between the minisuperspace vari- 
ables and their conjugated momenta. The authors found 
that noncommutativity modifies the cosmic dynamics at 
late times. Their result relies, however, on the study 
of an exact solution of the field equations that, in the 
commutative case, results in decelerating pace of expan- 
sion both at early and late times. It is clear that the 
latter is just a particular solution of the cosmological 
equations that can be attained only after a very care- 
ful arrangement of the initial conditions. Actually, as it 
was clearly and strictly shown in reference [23| - through 
the use of the standard tools of the (linear) dynamical 
systems analysis - the generic cosmological evolution at 
late times, in the case of a dilaton exponential potential- 
driven dynamics, depending on the region in the parame- 
ter space, can be governed either by an inflationary scalar 
field (SF)-dominated attractor solution, or by a scaling 
(also attractor) solution where neither the SF nor the 
background fluid dominate. In the latter case the ratio 
of the energy densities of the components of the cosmic 
mixture is a constant. As in the former case, at late times 
the expansion may be inflating even in this case. 

The results of the study in Ref . [22j , although relevant 
in what respects the question put forward: can noncom- 
mutativity affect the dynamics of the universe in the large 
cosmological scales?, are of limited importance in what 
regards their reach. Actually, their results can not be 
generic since these are based on a particular exact solu- 
tion. 3 In the present paper we aim at looking for a generic 
answer by the use of the standard tools of the dynamical 
systems analysis in the way the authors of Ref. |23j did. 
The basic steps of the formalism we use to build noncom- 
mutativity into our cosmological setting - as we already 
said it is based on the approach of reference [3| - are the 
following: i) starting from the Einstein-Hilbert action 
(plus a self- interacting scalar field action), the Hamilto- 
nian is derived for a Friedmann- Robertson- Walker space- 
time with flat spatial sections, ii) canonical quantization 
recipe is applied: the minisuperspace variables are pro- 
moted to operators and the WDW equation is written 
in terms of these variables, iii) noncommutativity in the 
minisuperspace is achieved through the replacement of 
the standard product of functions by the Moyal star 
product, and iv) semi-classical cosmological equations 



are obtained by means of the Wentzel-Kramers-Brillouin 
(WKB) approximation, 4 applied to the (equivalent) mod- 
ified Hamilton- Jacobi equation [25[. After this we choose 
adequate phase space variables and write the modified 
cosmological equations in the form of an autonomous 
system of ordinary differential equations (ODE). Then 
we follow the steps of Ref. [23| ■ 

In the present investigation we do not restrict our 
study to exponential potential only, but include a most 
general situation when a cosh-like potential is also con- 
sidered. Our results will indicate - in as much strict 
that the noncommutative effects not 
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terms as in 

only affect he early-times dynamics, as expected due to 
their quantum nature, but also the cosmic dynamics at 
late times. 

Before going any further, we want to mention the main 
"drawbacks" of our approach: i) we will be consider- 
ing noncommutativity in minisuperspace rather than in 
spacetime itself, and, ii) following the approach of Ref. 
4] we will be missing noncommutativity among the mo- 
menta conjugated of the minisuperspace variables. The 
former drawback is not as worrying as the latter one, 
since, although we are not dealing with direct space- 
time noncommutativity, the kind of noncommutativity 
in minisuperspace we will be considering, is expected to 
be a derived consequence of the former [4Ul4fl . Although 
the latter drawback is of more concern, however, our ap- 
proach here can be considered as a first step towards a 
more complete picture where noncommutativity of the 
conjugated momenta is also taken into account. 

We want to make emphasis in the fact that, given the 
above "drawbacks", and the "freedom" to choose min- 
isuperspace coordinates in order to build noncommuta- 
tivity into quantum cosmology 5 - a theory yet lacking 
a strict formulation -, the cosmological model studied 
here can be considered, at most, as a useful toy model 
to qualitatively study possible cosmological consequences 
of noncommutative effects. In this sense, this is just a 
first attempt at the latter goal, and, consequently, in the 
present paper we shall be concerned mainly with the cos- 
mological asymptotics and no attention will be paid to 
observational testing, a subject that deserves an indepen- 
dent publication. 

The paper has been organized in the following man- 
ner. In the next section [II] the basic mathematics of the 
model, as well as of the approach undertaken here, are 



3 In other words, the given solution may be unstable and very 
unlikely to arise in a real physical context. 



4 Also known as Jeffreys- Wentzel-Kramers-Brillouin (JWKB) ap- 
proximation |24| . 

5 Recall that, even if one follows the more standard assumption 
of noncommutativity among spacetime coordinates, one would 
not obtain the same results assuming noncommutativity, e. g., 
among cartesian or spherical coordinates. 
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exposed. We pay special attention to the following as- 
pects: i) the canonical quantization procedure for minisu- 
perspace variables, ii) introduction of noncommutativity 
through the Moyal star product, and, hi) derivation of 
the (semi-)classical cosmological equations by means of 
the WKB approximation. The main section lnTl is devoted 
to the study of the asymptotic structure of the modified 
(semi-classical) cosmological model, through the use of 
the dynamical systems tools. The exponential and cosh 
type of potentials are studied separately. Figures with 
the relevant phase space portraits are supplied. The re- 
sults obtained are discussed in detail in section IIV1 In 
section [V] the most important results of the paper are 
summarized and their reach commented on. We have 
added an appendix (section IVI[) where we discuss a pos- 
sible generalization of our results to include the presence 
of cold dark matter (CDM) in general instead of a scalar 
field fluid, through the exploration of a simple toy model. 
Here wc use units in which k 2 = 8irG = 1. 



S g = — j / d 4 xdt £(a, a), 



(2) 



where we have introduced the definition of the effective 
Lagrangean: C(a,a) = — 3e 3a d 2 . 

If, in addition to the standard EH (pure gravity) action 
S g , one introduces also a self-interacting scalar field (SF) 
action: 



where V(ip) is the self-interaction potential, while 
(dip) 2 = g^dptpdvip, then, the total effective Lagrangean 
under the integral in the resulting action: 

Stot = -^2 J ' d 3 xdt C tot (a,a,<p,(p), 
can be written in the following form: 



II. SET UP AND BASIC EQUATIONS 

As already mentioned in the introduction, here we 
follow the approach put forward for the first time in 
reference [4]. This time, however, we consider a ho- 
mogeneous and isotropic Friedmann-Robertson- Walker 
(FRW) spacetime with flat spatial sections. As a source 
of the Einstein's equations we choose a self-interacting 
scalar field (p. We will consider self-interaction poten- 
tials of the exponential and cosh-type. 



A. Hamiltonian approach 



£tot — 



„3q 



-3d 2 + -^ 2 - V(<p) 



(4) 



The Euler-Lagrange equations yield to the second Fried- 
mann (also Raychaudhuri) equation, and to the Klein- 
Gordon equation, respectively: 



2d + 3d 2 = --tp 2 + V(p), (p + 3aip = -V v , (5) 

respectively. Here, and in what follows, V tV = dV/dip. 

The next step is to go onto the equivalent Hamiltonian 
formulation. To this purpose we introduce the canonical 
momenta: 



To start with, we write the - 
Einstein-Hilbert (EH) action: 6 



general relativity (GR) 



(1) 



where k 2 = 8irG. If the homogeneous and isotropic FRW 
metric - flat spatial sections - is considered, then we can 
parametrize the metric as: ds 2 = —dt 2 + e 2a ^ 5ijdx l dx : > , 
where a(t) is the time-dependent scale factor. Then, 
after integrating by parts, the above action can be re- 
written as, 



6e 3Q d, 7iv 



dip 



e 6a ip, (6) 



respectively. The standard relationship between the clas- 
sical Hamiltonian % and the effective Lagrangean Ctot'- 
H(a, -K a ,ip,TT v ) = 7r Q d + ir v ip - C to t{a, d, ip, tp), yields to 



-3a 



H 



12 



[- 7 : 2 a + 6n 2 +12e ba V(i P )] 



(7) 



B. Mini-superspace variables and noncommutative 
WDW equation 



6 For a pedagogical and readable explanation of the procedure we 
are about to expose see section 2 of Ref. |2(I . 



Here, following reference Q, noncommutativity is 
achieved by an appropriate deformation of the usual 
(commutative) algebra of minisuperspace variables a, ip: 
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[ot,tp] = [a',<p'] = — z6», (8) 

where the new minisuperspace variables a' and tp' are 
noncommutative coordinates on a new (minisuperspace) 
base manifold |26[ , and the constant parameter 9 encodes 
the noncommutative effects. 7 

We think it is appropriate to make a comment on the 
physical motivations of the latter ansatz. Our assump- 
tion of noncommutativity among the expansion factor a' 
and the scalar field tp' is based in the results of previ- 
ous studies (see, for instance, references 043)- In our 
model the parameter 9 is expected to be a kind of effective 
(resulting) measure of noncommutativity amongst fields 
that are expanded in terms of the parameters related with 
noncommutativity of usual spacetime coordinates. The 
ansatz in noncommutative quantum cosmology is also in- 
spired in noncommutative quantum mechanics, due to 
the fact that a (a') and tp (ip ! ) play the role of "coordi- 
nates" in the minisuperspace, and appear as such in the 
WDW equation to study quantum cosmology (27| . 

Let us to promote the canonical momenta in ([7]) to 
quantum operators (d a = d/da, d v = d/dtp): ir a — > 
-?r a = —id a , ittp y Tttp = —id v , i. e., 

n U = ^ [dl - 6d 2 v + \2e^V{p)] . (9) 
The minisuperspace WDW equation then reads: 

H$(a,tp)=Q, (10) 

where ^ is the wave-function of the universe. 

Here we shall introduce the noncommutative effects 
through the Moyal star product of functions: 

U^(a, tp) = ->■ H * *(a, tp) = 0. (11) 

Since under the star operation the terms containing d 2 
and d 2 are unchanged, the effects of the Moyal/star prod- 
uct are reflected in the WDW equation through a shift 
in the potential: 

9 9 

V(a, tp) * *(a, tp) = V(a + -ft v , ip - -7r a )*(a, tp). (12) 

It has been shown that this is equivalent to performing 
the following change of variables |26j (see also [28j): 



7 In what follows, for definiteness, we consider 9 to be a non- 
negative quantity. 



a' = a + -Tt v , tp' = tp - -7t a , (13) 

where the involved field variables obey the corresponding 
commutation relationships in Eq. ((SJ). 

To make the point clear: we can write the WDW equa- 
tion modified through the replacement of the standard 
product of functions and operators by the Moyal/star 
product T-L * ^(a, tp) = 0, which results in Eq. (TUT) (sec 
below), or, equivalently, T-L^{a\tp') — 0, formulated in 
terms of the noncommutative coordinates a! , tp', of the 
new base manifold. Here we adhere to the former point 
of view, i. e., we keep the "old" commutative minisuper- 
space variables a, and tp, so that the noncommutativity 
will be encoded in the new terms containing the param- 
eter 9. The noncommutative counterpart of the WDW 
equation then reads: 

Wfltf (a, tp) = 0, =>• [d 2 a - %d% + U e (a, tp)}^> = 0, 

U e (a, p) = 12 e (e a -mo/a v ) v{ip + .6B {u) 

2 ooi 

If we further Taylor expand in the small parameter 
9 - this will entail assuming small 90S (a) /da and 
9dS(tp) / dtp in the equivalent Hamilton- Jacobi formula- 
tion (see below) - then, keeping up to linear terms: 

e C«"-3wa/a»0 = e 6 « (l - m A + o{9 2 )^ , 

^ + f^> = V'M + f^ + 0( ^ 

Recalling that we demand the operators d/da and 
d/dtp to act exclusively on the wave function "J, it can 
be shown that, in the above considered linear approxi- 
mation, the potential energy operator in equation (|14p. 
Ug can be written as follows: 

U e (a, p) = 12 e 6a (v - 3i9Vd v + jV,^d a ) . (15) 



C. WKB approach 

To get a consistent description of the cosmological dy- 
namics at any time we need a (semi) classical analog of 
the deformed WDW equation (|T4| . (p~5|) . Such a classical 
approximation to the obtained (modified) WDW equa- 
tion may be based on the Hamilton- Jacobi formulation 
and the WKB approach [U, [2||. Hence, following the 
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lines of Ref. [25| and references therein, lets assume sep- 
arability of the wave-function in its arguments a, tp: 

together with the following conditions on the derivatives 
of the S- functions: 



d 2 S(a) 



da 2 



fdS(c 



V dc 



dtp 2 



« m 



then, by substituting back into (O, (fT5|) . one obtains 
(recall that in the linear approximation considered here, 
one keeps terms up to 0(9j): 



(dS{a 



\ da 



dS(<p) 



12 e ba V + 39V- 



dtp 



2 V >* 



dS(a) 
da 



.(16) 



Next, in order to get the classical equations, as it is cus- 
tomary in the standard WKB procedure, we make the 
following identifications: 



D. Cosmological equations 

Equations (fI5]l and (I2TJ1) form a closed system of differ- 
ential equations. The Raychaudhuri (also, second Fried- 
mann) equation can be derived by taking the time deriva- 
tive of the Friedmann equation (I18|) and substituting tp 
from the modified Klein-Gordon equation (I2U1) . We ob- 
tain 



3d 2 = --p 2 + V + 39e 3a (aV ip + pV). (21) 



2a 



Equations (jig]). (|20|). and (|2l|). are the modified GR equa- 
tions of motion provided by the WKB approach, for the 
case where 9 is a small quantity. 

A straightforward inspection of the modified cosmo- 
logical equations show that the effect of the kind of non- 
commutativity explored in this paper is to modify the 
effective self-interaction potential of the scalar field. Ac- 
tually, the second and third terms in the RHS of Eq. (f2Tj) 
can be joined together into the effective potential V e ff- 



V eff = V + 36 e 3a {dV v + <pV), 



(22) 



dS(a) 
da 



dS(<p) 
dtp 



(17) 



where 7r Q and tt v , are just the canonical momenta given 
by Eq.©. Hence, after substituting back (IT7l) into Eq. 
(fT6|) (taking into account the definitions of the canonical 
momenta in (Eq. ©), one obtains the following modified 
Friedmann constraint: 



3d 2 = -ip 2 + V + 39 e 3a (aV v + tpV) . (18) 

A dynamical equation to determine tp can be derived 
coming back to the Hamiltonian constraints (|14l) . (IT51) . 
which, in terms of the classical momenta 7r Q , n v , result 
in: 



which means that we can define the effective parametric 
energy density and pressure of the scalar fluid: p e J' = 
tp 2 /2 + V e ff, p e Jf = p 2 /2 — V e ff, respectively. Hence, 
the following set of " standard" field equations for a SF- 
sourced cosmology is obtained: 



3d' 



= p°ff 



2d 



P 



eff + 3a{p% sl 



f 3d 2 = 
P e J f ) = 



-P e / f , 
0. 



An alternative interpretation is possible however. Ac- 
tually, it is really tempting to identify the following effec- 
tive energy density and pressure, respectively, of a "non- 
commutative" (NC) vacuum fluid (see Ref. [29( for a 
related model): 



= - P e = 39 e Aa (aV v + pV) 



(23) 



n e = 



12 



-(-7T 2 +67T 2 ) + 



„3q 



V + 36Vw v --V, l/> v a ). (19) 



Here the noncommutative effect is encoded in the small 
parameter 9. Then, by taking into account Eq. ([6]), the 
canonical Hamilton's equation, -k v — —&Hg/dtp, yields 
to the desired equation: 



ip + 3atp = -V v - 39e 3a (aV iVV + tpV v ). 



(20) 



After the above identification one may re-write the cos- 
mological equations (TTB")) . ([2U|) . (|2"Tj) . in the following com- 
pact form (we also write the continuity equation for the 
effective noncommutative fluid): 



3d = p v + p g , 2a 
p v + 3a(p ip +p<p) = 



3d 2 = -p v - pg, 
-<PPe, Pe = <PPe, 



(24) 



where, as usual, p v — tp 2 /2 + V, is the (parametric) 
energy density of the scalar field fluid, while p v — tp 2 /2 — 
V, is its parametric pressure. 
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Notice that, although the conservation equation is not 
obeyed by each separate component of the cosmic mix- 
ture, the mixture, as a whole, actually does obey the 
continuity equation: p to t + 3d (/?t t + ptot) = 0, where 
Ptot — Pip + Pe, and p to t — Pip + Pe- It is evident that, 
written in the above suggestive form, the cosmic dynam- 
ics is governed by additional, non-gravitational interac- 
tions between the scalar field fluid and the effective NC 
fluid, through the source term tppg. The dynamics of in- 
teracting cosmic fluids is a very interesting scenario to 
look for solutions to the cosmic coincidence problem (30l — 

Another interesting feature of this alternative interpre- 
tation is that the effective NC fluid behaves just like a 
cosmological "constant" (pe = —pe), where, due to the 
source term ippe in the continuity equations, the cosmo- 
logical constant is in fact a dynamical quantity. This 
means that the impact of noncommutativity on the cos- 
mic expansion at late times may be non- vanishing. 

In the next sections we shall explore the asymptotic 
properties of the above cosmological model, to look for 
generic features of the impact noncommutativity (of the 
kind considered here) might have on the large-scale cos- 
mic dynamics. 



III. DYNAMICAL SYSTEMS 

Our task here will be to write the equations (fT8|) . (|20|l . 
H]), in the form of an autonomous system of ordinary 
differential equations. To this end we have to choose 
appropriate phase space variables. In the present case 
our starting choice is the following: 



-3x(l-ir 2 ) + 

y' = - x 2 ), z = 



f + 3y(y/6x - zT) 
1 + 3y(V6x — z) 

V6xz 2 (T- 1), (27) 



where T = VV^/V v , and the tilde denotes derivative 
with respect to the time-ordering variable a = In a (ba- 
sically the number of e- foldings) . 

The above system of ODE is a closed system of ordi- 
nary differential equations only if the parameter Y can 
be written as a function of the variable z\ V = T(z). In- 
deed, for a large class of self-interaction potentials which 
are very popular in the cosmological context, T can be 
written in the form of a polynomial in z. 

The phase space relevant for this problem depends on 
the particular kind of potential considered, but, in gen- 
eral it can be defined as: 



$ = {(x,y,z)\y>0, z>0...}. 



(28) 



The dots mean additional constraints on the phase space 
variables, coming from the particular functional form of 
the self-interaction potential under consideration. We 
choose y > since we are interested in expanding cos- 
mologies only (d > 0, recall that we have chosen 9 to be 
a non- negative quantity). 

Several magnitudes of cosmological relevance are the 
following. The deceleration parameter (q = — (1+d/d 2 )): 



s/Ec 



y = 8e 3a a, 



V ~f- (25) 



Notice that x amounts to the dimensionlcss kinetic en- 
ergy of the scalar field squared, while y encodes the non- 
commutative effects. The phase space variable z is sen- 
sitive to the kind of self-interaction potential chosen. In 
what follows we restrict our study to non-negative z-s, 
since the negative sector can be obtained from the non- 
negative one by the transformation ip — > — cp. 

The Friedmann equation (jTHJ) yields to the following 
constraint relating the dimensionless SF potential energy 
with the variables x, y, z: 



V 

3d2 



1 - X 1 



1 + 3y{V&x - z) ' 



(26) 



The autonomous system of ODE that can be obtained 
out of the cosmological equations (TT5j) . (|2"0)) , (f2"Tj) . is the 
following: 



-1 + 'ix 2 



(29) 



and the dimensionless energy densities of the NC fluid 
and of the SF: 



Pe _ 3(1 - x 2 )(V6x - z)y 
3d 2 ~ 1 + 3y(V6x - z) 



, n v = l-Qg, (30) 



respectively. Recall that, since the effective NC fluid be- 
haves like vacuum (jpe = —pe), then the effective equation 
of state (EOS) parameter, ujq = —1, always. 

In general, for potentials which vanish at the minimum 
- usually associated with relevant late-time behavior - 
the phase space variable z is undefined, and one have to 
choose an alternative variable. One possible choice is the 
following: 



(31) 



z+l 
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where, since z is a non- negative quantity, then < z < 
1. Besides, the resulting phase space spanned by the 
variables x, y, z is noncompact only in the y-direction, 
so that several critical points might scape our study. A 
way out is to transform also the variable y: 



1 



1 



y 



l 



y 



y 



y 



After the above transformations the phase space spanned 
by x, y, z, is a compact space, which means that all of 
the existing equilibrium points are confined to a bounded 
region. 

In terms of x, y, z, the autonomous system of ODE 
(l27l) transforms into: 



l-z 



-3z(l-x 2 ) + A /~(l-z 2 ) 



x [ yz + 3(1 - y) [VExz - (1 - z)T(z)] 

y' = -3y(l-y)(l-x 2 ), 
z = %/6x(l - zf[T(z) - 1]. 



(32) 



As already said, the phase space is now a compact 3- 
dimensional space spanned by the coordinates x, y, z, 
and can be defined in the following way: 

$ = {(x,y,z)\-l < x< 1, 0< y < 1, 0< z < 1}. (33) 

Recall that the ranges y > 1 and y < 0, are not be- 
ing considered since we are interested in describing the 
dynamics of expansion (contracting phases are not of in- 
terest to us). 

Such relevant cosmological parameters as the SF and 
NC dimensionless energy density parameters, together 
with the scalar field EOS parameter, can be given in 
terms of the new variables as follows: 



and 



3(1 - x 2 ){l - y)[(VGx + l)z - 1} 
yz + 3(l-y)[(V6x + l)z-l] ' 



£l v = l 



1 



2(x 2 - l)yz 



respectively, 
changed: q = 



yz + 3x 2 (l - y)[(V6x + l)z - 1] ' 
Other relevant magnitudes remain un- 

-l + 3x 2 , LUg = -1. 



A. Exponential potential 

For the particular case when z — A =>- z 
is a constant: V(<p) — Voexp(— Xip) 



= V(A + 1), 
r = 1. In 






FIG. 1: Phase portraits for the exponential potential case 
(subsection MI All for different choices of the exponential pa- 
rameter A: i) A — (top left), ii) A = 0.5 (top right), and iii) 
A = 3 (bottom). In this, as well as in the subsequent figures, 
the x, y (and z) axes span the compact phase space variables 
x, y (and z) respectively. Several orbits corresponding to dif- 
ferent sets of initial conditions are shown. In each case the 
initial conditions on the "noncommutative" coordinate y cor- 
respond to either j/(0) = 0.75 (upper groups of orbits), or to 
y(0) — 0.25 (lower groups of orbits) only. The attractor struc- 
ture of the scalar field/noncommutative fluid scaling solution 
(critical point P v /e) when A 2 < 6 (upper panels), is evident. 



this case the system of ODE (|3"2"|) appreciably simplifies, 
and instead of a 3-dimensional phase space, one have a 
2-dimensional one, spanned by the variables x, y: 



x' = -^(l-x 2 )(V6x-X), 

y' = -3y(l-y)(l-x 2 ). (34) 

We have, also, that q = — 1 + 3x 2 (as before), and 

3(1 - x 2 )(V&x - - y) 



VLB = 



y + 3(VSx-\)(l-y) 



The main qualitative properties of the phase space 
spanned by the variables x, y, z, for this case, are sum- 
marized below. 
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1. Critical points: commutative (GR) limit 

In the commutative GR limit (y = =>■ y = 1), 
the phase space is a linear space, and the critical points 
are x = ±1 and x = X/\&. Perturbing (linearly) in the 
neighborhood of these points one obtains: 

respectively, where e is the small perturbation. Hence, 
the perturbation evolves according to e±(a) = eo exp((6=F 
y/6X)a), for the points x — ±1, and e(a) = eoexp((A 2 — 
6)a/2), for x — A/v6) respectively. This means that, 
while the SF kinetic energy-dominated solution corre- 
sponding to the choice x = 1 (also stiff- fluid solution), 
can be stable for A > y6j the one corresponding to 
x = — 1 is always unstable. For A 2 < 6 the SF ki- 
netic/potential energy-scaling solution x = A/V6 is sta- 
ble, otherwise it is also an unstable critical point. 

2. Critical points: general case (includes 
noncommutativity ) 

In the general case (including noncommutativity), the 
critical points, Pj = (x,i,y~i), of the above autonomous 
system of ODE are: 

• SF kinetic energy-dominated solution: P£ — 
(±l,j/). The cosmological parameters, in this case, 
are: 

ile = 0, = 1, q = 2. 

Since, in both cases, there is one vanishing eigen- 
value of the corresponding linearization matrices: 
i) Ai = 0, A 2 = 6 - V6A, for P+, ii) Ai = 0, 
A2 = 6 + \/6A, for P^, these critical points are 
non- hyperbolic . 

• SF scaling solution: 

Pk/v = (A/VM). The di- 
mensionless density parameters, and the decelera- 
tion parameter are equal to: 

A 2 

n g = 0, n v = 1, q = -1+ y. 

The eigenvalues of the linearization around this 
point are: Ai = 3 - A 2 /2, A 2 = -3 + A 2 /2, so 
that it is a saddle equilibrium point . 

• SF/NC fluid-scaling solution: P v/e = (A/-\/6,0). 
The relevant cosmological parameters are: 

A 2 A 2 A 2 

n e = i-— , n v = — , g = -i + _ , 

6 6 2 



while the eigenvalues of the linearization matrix 
are: Ai = A 2 = -3 + A 2 /6. For A 2 < 6 this is 
the future attractor in the phase space. 

Summarizing the results of this subsection: i) the non- 
commutative effects modify the early-times cosmic dy- 
namics by erasing any past attractor (see the figure [lj , 
ii) these effects also modify the late-times dynamics by 
replacing the two possible future attractors according to 
the standard GR-limit (either the stiff-fluid solution, or 
the SF kinetic/potential energy-scaling solution) by the 
SF/NC-scaling solution (critical point P v /e)- 

B. Potentials of the cosh-like type 

Here one is considering potentials of the following kind: 

V(ip) = V Q [cosh(A^) - if , (35) 

so that 

pXsmh(Xip) rfl = , 1 pX?_ 
cosh(A(^) - 1 ' K ' 2p 2z 2 ' 

Hence, in terms of the variable z: 

_ (2p-l)(l-z) 2 +p 2 A 2 z 2 
[Z) ~ 2p(l - zf 

besides; 

(p 2 A 2 -l)z 2 + 2z-l 

r(z) - 1= — w^w — • 

If we substitute back into (|32|) , and then we look for the 
equilibrium points of the resulting autonomous system of 
ODE, we obtain the results which we list below (see also 
the tables HUH andHU). 



1. Critical points: commutative (GR) limit 

The commutative case corresponds to the choice y = 1, 
so that y = => 9 = (the case a = is of no interest 
for cosmology). The corresponding critical points and 
their properties are summarized in table HI Here we list 
these results. Recall that critical points are given by pairs 
(x,z). 

• Kinetic energy-dominated critical points: 
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TABLE I: Critical points and their properties for the cosh-like potential V(ip) = Vo [cosh(Ay) — l] p : the commutative GR limit. 
The eigenvalues of the corresponding linearization matrices are shown in the right-hand columns. 

C. Points x z Existence ui v q Ai A2 

1 Tfp\ Always 1 2 6 - V6pX \/6A 

Pjc — 1 TTpA " 1 2 6 + ^pA -y^A 

P P x 1 » 1 1 V 2 >? 1 1 V 2 >- 2 o 1 p 2 A 2 \2 



-1 1 -§(l + ,/l-§ M 3) -I l-Jl-§^ 



TABLE II: The critical points of the autonomous system (|32[) for the cosh-like potential V(<p) = Vo [cosh(A(^) — general 
case (includes noncommutativity). Here we have defined r = (2p — 1 + y / p 2 A 2 (l — 2p))/ (2p — 1 +p 2 A 2 ). 

C. Points x y z Existence fi w fie w B g 

Pt 
Pk 

Pk/v 



Pnc 



1 


V 


1 

1+pA 


Always 


1 





1 


2 


-1 


.'7 


1 

1+pA 


?) 


1 





1 


2 


J-'A 
\/6 


1 


1 

1+pA 


)? 


1 





2 \ 2 

_1 + £_*_ 
1 T 3 


2 \ 2 

1 T 2 





1 


1 


?) 


1 





-1 


-1 





1 





n 


Undefined 


Undefined 


-1 


-1 








r 


p< 1/2 





1 


-1 


-1 



also known as stiff-fluid solution. For these points 
ui v = 1, and q = 2, while the eigenvalues of the 
linearization matrix are Ai = 6 — \/6pA, and A2 = 
v6A, for Pi, while for PZ these are: Ai = 6 + 
\/6pA, A 2 = — v6A. Hence, for pX < v6 ; fjf is 
the past attractor in the phase space, while, for 
pX > \/6, it is a saddle critical point instead. The 
critical point P^ is always a saddle point in the 
phase space. 

• SF kinetic/potential energy-scaling solution: 

Pk/v = (^TTp->)- 

For this case the SF EOS and the deceleration pa- 
rameters are given by cj^ = —1 + p 2 A 2 /3, q = — 1 + 
p 2 A 2 /2, respectively. Since the eigenvalues of the 
linearization matrix for this case Ai = — 3+p 2 A 2 /2, 
A2 = pX 2 , then this critical point is a source (past 
attractor) whenever p 2 A 2 > 6, while for p 2 A 2 < 6 
it is a saddle equilibrium point. For p 2 A 2 < 2 the 
corresponding cosmological solution describes ac- 
celerated expansion. 



• SF potential energy (V)-dominated solution: Py = 
(0, 1). The scalar field fluid mimics a cosmological 
constant since uj v — — 1 and the peace of the cosmic 
expansion is accelerated (q = —1). This solution is 
always the future attractor since the real parts of 
the eigenvalues of the linea rization mat rix are both 
negative: 2Ai, 2 = -3(1 ±^1 - 2pA 2 /3). For pX 2 > 
3/2 the future attractor is a spiral equilibrium point 
(see the Fig. EJ. 



Summing up: within the framework of standard (com- 
mutative) GR-limit, the past attractor can be either the 
stiff- fluid solution if pX < v6 (the K/V energy-scaling 
solution and the conjugated stiff fluid solution are both 
saddle critical points), or the K/V energy-dominated so- 
lution if pX > i/6. In the latter case the stiff fluid solu- 
tion (and its indistinguishable conjugated, in this case) 
is always a saddle equilibrium point. The future (late- 
time) attractor - a spiral point for pX 2 > 3/2 - is always 
the inflationary SF potential energy-dominated solution. 
These features are illustrated in the phase portrait in the 
figure [H 
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TABLE HI: Eigenvalues of the line arization matrices corresponding to the critical points in TabHTl Here we have defined: 
m = 3y/l - (VI - 2p-p 2 A 2 ) 2 /6p 2 . 



C. Points x y z Ai 



P+ 1 y 6-V6pA V6A 

P~ -! y -L- o 6 + V6pA -V6A 

Pv Oil 3 -3 + m -3-m 

Py 10 3 -3 + m -3-m 

Pnc r -3 -3 + - 3pA 2 -3 - ^9 - 3pA 2 



2. Critical points: general case 

While in the GR-limit there are found four critical 
points, in the general case - considering NC effects of 
the kind explored here - there are found two addi- 
tional critical points. Although all these critical points, 
Pi = (xi,yi,Zi), are exposed in the table Ull (the eigenval- 
ues of the corresponding linearization matrices are shown 
in the table here we list them and discuss their basic 
properties. 

• Kinetic energy-dominated critical points (stiff fluid 
solution) : 

^ = ( ±1 ' 5 'TT^)- 

For these points fi^ = 1 (fig = 0), uj v — 1, and 
q = 2. For P^ the eigenvalues of the linearization 
matrix are Ai — 0, A2 = 6 — V^pX, and A3 = v6A, 
while, for P^ these are: Ai = 0, A 2 = 6+\/6pA, and 
A3 = — v6A. In both cases, since one of the eigen- 
values is vanishing, these are non-hyperbolic points 
and we miss part of the information we could ob- 
tain from the application of the standard tools of 
the (linear) dynamical systems analysis. However, 
for P£ , as long as pX > ^/6, this is a saddle equilib- 
rium point since A2 and A3 are of opposite sign. In 
the case P^ , since for positive p the eigenvalues A2 
and A3 are of opposite sign, then we can say with 
certainty that this is always a saddle critical point. 
Additional information can be extracted from the 
inspection of the phase portrait (see Fig. |3]). It can 
be corroborated that there are no past attractors 
there. 



• SF K/V energy-scaling solution: 

^ = (S' 1 'TT^)' 

The cosmological magnitudes of relevance are given 
by: fi v = 1 (fig = 0), u v = -1 +p 2 A 2 /3, g = -1 + 
p 2 X 2 /2. The eigenvalues of the linearization matrix 
for this case are: Ai = pX 2 , A2 = —3 + p 2 A 2 /2, 
A3 = 3 — p 2 X 2 /2. Since the eigenvalues A 2 and A3 
are always of opposite sign, then this critical point 
is always a saddle point in the phase space, unlike 
in the GR-limit, when this point can be also a past 
attractor. 

• SF V-dominated solution: P v = (0, 1, 1). The rel- 
evant parameters are Q v — 1 (fig = 0), uj v = — 1, 
and q = — 1, while the eigenvalues of the lineariza- 
tion matrix are: Ai = 3, A2,3 = — 3±m (the quan- 
tity m is defined in the heading of Tab. . This 
is always a saddle equilibrium point in the phase 
space. 

• Conjugated SF V-dominated critical point: Py = 
(0,1,0). In this case both fi v and fig are unde- 
fined. The other quantities remain the same as in 
the former case, including the eigenvalues of the 
linearization matrix. This point can be associated 
with the minimum of the cosh- like potential (I35[) . 

• NC-dominated solution: Pnc = (0, 0, r), where the 
quantity r has been defined in the heading of Tab. 
III1 This solution is dominated by the energy den- 
sity of the effective noncommutative fluid ilg = 1 
(fi„ = 0), and exist whenever p < 1/2. It is as- 
sociated with accelerated expansion since q = — 1. 
Whenever it exists, it is the late-time (future) at- 
tractor since the real parts of the eigenvalues of 
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FIG. 2: Standard GR-limit for the cosh-like potential. The 
following values of the free parameters have been chosen: 
A = 5, p — 0.25. Several orbits corresponding to different 
sets of initial conditions are drawn. The SF potential energy- 
dominated solution Py = (0, 1) is the late-time attractor. It 
is seen that the orbits of the autonomous system of ODE 
spiral down to the future attractor, which is associated with 
coherent (damped) oscillations of the cosmological scalar field 
around the minimum of the potential. 




the corresponding linearization matrix are nega- 
tive: Ai = -3, A 2 , 3 = -3±V9~3pA 2 . For p\ 2 > 3 
this is a spiral critical point since the eigenvalues 
of the linearization matrix are complex numbers. 

A brief summary of the main results can be given: i) 
noncommutative effects modify the early-times dynam- 
ics by erasing any past attractor (neither the stiff-fluid 
solution, nor the SF K/V-scaling solution can be source 
points), ii) noncommutative effects also modify the late- 
times dynamics: the stability of the SF potential energy- 
dominated solution (the future attractor in the standard 
GR-limit) is modified: it is now a saddle critical point. 
For p < 1/2 the future attractor in the phase space 
is the inflationary NC-dominated solution (equilibrium 
point Pnc)- Additionally, the way in which the orbits 
in the phase space approach to the late-time attractor is 
also modified: while in the standard GR-limit the future 
attractor is a spiral point for p\ 2 > 3/2, in the general 
case, due to the influence of the noncommutative effects 
of the kind considered here, this attractor is a spiral point 
for pX 2 > 3. The above enumerated features can be ap- 
preciated in the figure [3l 

IV. DISCUSSION 

The question that stands in the title of this paper is 
not as trivial as it seems. To start with, noncommuta- 
tivity, if it really have played any roll in the cosmolog- 



FIG. 3: General case for the cosh-like potential (includes non- 
commutative effects). The flux in time a of the correspond- 
ing autonomous system of ODE is shown for A = 5, p — 0.25 
(top). The projections of the phase space onto the different 
phase planes are also shown (bottom) . While the existence of 
the future attractor - critical point Pnc = (0, 0, 0.36) - corre- 
sponding to the solution dominated by the noncommutative 
effects (see the table (TTJ , is evident, there are not found past 
attractors. Due to our choice of the free parameters above, 
pA 2 > 3, the late-time attractor is a spiral critical point. For 
p > 1/2 no future attractor can be found in the phase space 
neither. 



ical dynamics, is expected to have influenced the very 
early stages of the cosmic expansion, when, presumably, 
quantum effects had an appreciable impact on the grav- 
itational interactions of matter. However, the relatively 
recent discovery that our universe expands at an accel- 
erated peace, has put forward the possibility that tiny 
quantum effects that might have passed unnoticed in the 
recent past of the cosmic history, might be the cause of 
the present speed up of the expansion. Actually, the sim- 
plest model that accommodates the recent acceleration 
of the cosmic expansion rests on the anti-gravitating ef- 
fect of the quantum vacuum (the cosmological constant). 
Since the energy density of the vacuum (p vac oc A) does 
not evolve with the expansion, then, even if it has been 
a tiny fraction of the matter-energy content of the uni- 
verse in the past, as long as the remaining components 
of the cosmic mixture dilute with the expansion, there is 
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a moment in the future (present) of the cosmic history 
when the negative pressure of the vacuum starts domi- 
nating the cosmic dynamics, resulting in a new period of 
inflation (also referred to as late-time inflation). 

Therefore, it makes sense to ask whether noncommu- 
tative effects can have any appreciable impact on the des- 
tiny of the cosmological evolution either. Besides, there is 
no evidence that its influence on the early-times dynam- 
ics is a generic feature. This is why we have focused in 
the study of the asymptotic properties of FRW cosmolog- 
ical models in connection with noncommutative quantum 
cosmology. 

Here we have concentrated in a semi-classical WKB ap- 
proximation to the minisuperspace WDW equation gen- 
eralized to encompass noncommutativity of field vari- 
ables a and ip by the introduction of the Moyal star prod- 
uct. The resulting FRW cosmological equations ([18]). 
(|20|) . (|2Tj) can be given an (attractive) alternative inter- 
pretation in the form of the equations ([M]) . Written in 
the latter form, the noncommutative effects of the kind 
considered in this paper, may be encoded in an additional 
(effective) NC matter term which equation of state tracks 
that of the cosmological constant pg = —pe, so that, 
its possible inflationary effect is envisioned. However, 
unlike the usual cosmological constant term, the energy 
density of the NC effective fluid does actually evolve as 
the cosmic expansion proceeds, thanks to additional non- 
gravitational interactions with the scalar field through a 
source term oc <ppe- s These interactions are switched on 
by the Moyal star product, which mixes the metric and 
the scalar field components. 

Since, in general, the equations ([18]). (|20|) . (f2Tj) (or 
their equivalent (|24|)) are very difficult to solve analyti- 
cally, an alternative way around is to invoke the dynami- 
cal systems tools to extract very useful information about 
the asymptotic properties of the model instead. In this 
regard, knowledge of the equilibrium points in the phase 
space - corresponding to a given cosmological model - is 
a very important information since, independent on the 
initial conditions chosen, the orbits of the corresponding 
autonomous system of ODE will always evolve for some 
time in the neighborhood of these points. Besides, if the 
point were a stable attractor, independent of the initial 
conditions, the orbits will always be attracted towards 
it (either into the past or into the future). Going back 
to the original cosmological model, the existence of the 
equilibrium points can be correlated with generic cosmo- 



logical solutions that might really decide the fate and / or 
the origin of the cosmic evolution. 

What the results of our dynamical systems analysis 
have revealed is that, independently of the kind of self- 
interaction potential considered: i) exponential potential, 
or ii) cosh-like potential, the noncommutative effects of 
the kind considered here 9 affect not only the early-times 
dynamics but also modify the late-time behavior. Below 
we will discuss this issue in detail. 



A. Exponential potential 

Exponential potentials and their combination have 
been intensively studied in the recent past in connec- 
tion with cosmological applications [33j]. These arise in 
supergravity and in supcrstring after dimensional reduc- 
tion |34l |. 

The dynamical systems study of the NC cosmological 
model considered here reveals that, while for standard 
(GR) FRW cosmology the stiff fluid solution, x = 1, 
is always the past attractor and the scaling solution, 
x = A/V6 =>■ ip 2 /2V = A 2 /(6 - A 2 ), is the future 
attractor, 10 in the most general case, which includes non- 
commutative effects, there are no past attractors at all. 
In other words, the noncommutative effects modify the 
early-time dynamics by erasing any past attractor from 
the phase space. This means, in turn, that the starting 
point of the cosmic dynamics is uncertain: if we evolve 
the cosmological equations from given initial data (given, 
say, in the present epoch) back into the past, the re- 
sult is highly dependent on these data. However, the 
stiff fluid solution continues being an equilibrium point, 
which means that the cosmological evolution may evolve 
for some time in the vicinity of this solution. Hence, 
the noncommutative effects modify the stability proper- 
ties of the equilibrium configuration associated with the 
early-times cosmic dynamics. 

The surprising fact was to find that the NC effects also 
modify the late-times dynamics. Actually, in the gen- 
eral case, the SF kinetic/potential energy-scaling solu- 
tion, P K/V = (A/%/6, 1) 6 = 0, <p 2 /2V = A 2 /(6-A 2 ), 
is always a saddle critical point (it was the late-time at- 
tractor in the GR- limit), while the late-time attractor is 
the SF/NC fluid-scaling solution P v/g = (A/v^O) 



Solar system measurements impose severe constraints to occur- 
rence of additional non-gravitational interactions of matter, how- 
ever, when these interactions involve dark components of the cos- 
mic mixture, as in the present case, the issue is subtle and, in 
general, these constraints may be evaded 31.. 1 32H . 



We frequently repeat this sentence because there are several dif- 
ferent ways to build noncommutativity into a given cosmological 
setting. 

The case for A > \/6, where the stiff fluid and the scaling solu- 
tions exchange their stability properties, is not being considered 
since this condition would imply negative energy. 
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Q, v /VLq = A 2 /(6 — A 2 ). This solution is inflationary (i. 
c, it is correlated with accelerated expansion of the uni- 
verse) whenever A 2 < 2. 



B. Cosh-like potential 

The cosh-like potential ([331) . 

V(ip) = V [cosh(A^) - if , 

has been studied in connection with the so called scalar 
field dark matter (SFDM) models 

In a natural scenario for cosmic dynamics, the scalar 
field tp runs from arbitrarily large negative values (\ip\ 3> 
1) to vanishing ones (\ip\ <C 1)- In consequence at early 
times the dynamics is driven by an exponential potential 



p A /2), ii) the SF potential energy-dominated solution 



1 



Vo 



M » J V(<p) « -f e-^, 

whereas at late times it is associated with a power-law 
potential: 



M « ~ V(<p) « im 2 ^, 



^oA 2 



In general, for positive p-s, there is an oscillatory phase 
around the minimum of this potential which plays an 
important role in the late-time dynamics. It has been 
demonstrated [38| that in the limit when the oscillation 
period is much smaller than the time scale of the cos- 
mic expansion, coherent scalar field (damped) oscillations 
behave like a fluid with (p v ) — (uj^) p lp) where the mean 
equation of state of the fluid (u> v ) depends upon the form 
of the scalar field potential V(ip) [38|. In particular, for 
V = V (p n (uj v ) = (n-2)/(rc + 2). Hence, for the po- 
tential (|55|) . at late times - during the oscillatory phase 
- the mean equation of state is determined by 



Pv 



p-1 



P 



1 



Notice that for p = 1 the scalar field behaves like pres- 
sureless dust, (ljJ) = 0. A scalar field potential with this 
value of p could therefore play the role of cold dark matter 
(also known as SFCDM) in the universe. For, p < 1/2, 
this potential is a good candidate for quintessence models 
of dark energy (DE) jaEH- 

The known results of the dynamical systems study of 
this kind of potential - within standard, commutative, 
FRW cosmology - show that (see subsection MI B ip : i) 
The past attractor can be either the stiff fluid-dominated 
solution P£ = (1,1/1 +p\) {u) v = 1, q = 2), if pX < 
\/6, or, whenever, p\ > ^/E, it is the SF scaling solution 
Pk/v = (pA/>/6, 1/1+pA) (u v = -l+p 2 A 2 /3, q = -1 + 



Pv = (0, 1) (w ¥ 



T, meaning that ip = tpo, 



V((f) — Vq), is always the past attractor. For pX 2 > 3/2 
this is a spiral point signaling that the field ip performs 
damped (coherent) oscillations around the minimum of 
the potential V(<p) at late times, until the stable de Sitter 
solution is attained. This oscillatory stage is what can be 
identified properly with SFCDM [3jj . The above features 
are clearly illustrated in the figure [5J 

What kind of modifications of the above picture are 
produced by the noncommutative effects? According to 
the results of subsection MI B 21 after NC effects of the 
kind explored here - WKB approach to the minisuper- 
space WDW equation, supplemented with the Moyal star 
product - are switched on, the past asymptotic structure 
of the phase space is modified: there are no past attrac- 
tors there. 

The interesting finding is that the future asymptotics 
of the phase space is also modified. Actually, first, the 
stability of the SF potential energy-dominated solution 
(critical point Py above in subsection MI B 2]) . which was 
the late-time attractor in the standard commutative the- 
ory, is modified: it is now a saddle equilibrium point. 
Second, an additional conjugated (twin) SF potential 
energy-dominated solution P v arises. Third, the oscilla- 
tory behavior of the perturbation is also modified by non- 
commutativity: while in the standard GR picture coher- 
ent oscillations arise for pX 2 > 3/2, in the NC-modificd 
picture the damped oscillations occur whenever pX 2 > 3. 
Finally, the late-time attractor is the inflationary solu- 
tion that is dominated by the NC effects (critical point 
Pnc for which Qg = 1, see Fig. |3J). The latter attractor 
exists whenever p < 1/2. Hence, the noncommutative 
effects may be the cause that the universe inflates at late 
times. 

A genuine objection against our finding that the 
present speed up of the cosmic expansion might be due to 
the effect of noncommutativity, can be based on the argu- 
ment that the scalar field itself (alone) may fuel the late- 
time inflation. So that, why to complicate that simple 
quintessential picture with the inclusion of unknown non- 
commutative effects of quantum nature originated very 
early in the cosmic history? 

A reply to this kind of objections can be based on 
the following arguments. First, the SF models of dark 
energy are plagued with serious problems: fine tunning of 
initial conditions, the coincidence problem, etc. Besides 
the notion of dark energy itself is very unappealing and 
faces serious challenges to find support in the standard 
scheme of the fundamental interactions. 11 



This is why alternative ways to explain the present stage of ac- 
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A second argument is a bit more technical. Recall 
that a scalar field with a cosh-like potential of the kind 
(|55|) , can be a nice model of dark matter if one chooses 
p = 1 [35l - [37j . while, for p < 1/2, it is a good candidate 
for quintessence model of dark energy instead [H, |39| . 
However, the latter statements are true only for non- 
interacting scalar fields - no additional interaction with 
the other components of the cosmic mixture. In par- 
ticular, in the present case where the scalar field inter- 
acts with the effective NC fluid, the above statements 
are not true. This can be clearly seen by recalling that 
the influence of noncommutativity can be alternatively 
understood as a modification to the SF self-interaction 
potential. Hence the range of values of the parameter 
p for which the cosh-like potential correctly explains the 
dark matter, is shifted to lower values. Even if under 
the effects of noncommutativity, the (interesting for cos- 
mology) picture with a late-time attractor inflationary 
solution, arises for p < 1/2, due to the mentioned shift 
in p, the scalar field with the potential (|35|) can be a 
good candidate to account for the dark matter. This 
statement is supported by the results of the dynamical 
systems study discussed above: for, p\ 2 > 3, coherent 
(damped) oscillations of the SF perturbation around the 
NC-dominated solution, arise. This oscillatory behavior 
is what can be interpreted as the (SF)CDM [33]. Hence, 
to close our line of argument; in our model the scalar field 
plays the role of CDM, while the noncommutative effects 
account for the late-time speed up of the expansion of 
the universe. 

In the Appendix, through the inspection of a simplified 
toy model, we avoid relying on the scalar field component 
to show that the noncommutative effects can, in fact, be 
a nice candidate to explain the cosmic speed up at late 
times. 



V. CONCLUSION 

In the present paper we have addressed the question, 
formerly explored in Ref. [22[ , about the possible impact 
of noncommutative effects of quantum nature on the dy- 
namics at large cosmological scales, by the study of a 
simplified model 0, [H, Ell • 

While in 22] the issue was investigated in connec- 
tion with a given particular solution of the corresponding 



celeration of the cosmic expansion have been explored. These 
include modifications of the laws of gravity. It is in this latter 
vein where our proposal fits. Modifications of gravity of quantum 
nature, such as inclusion of noncommutative effects, can be an 
interesting alternative to explain this — up to date — mysterious 
speed up of the expansion of the universe. 



modified cosmological equations (and for an exponential 
dilaton/scalar field self-interaction potential), here we 
have approached this subject from the dynamical sys- 
tems perspective, and, additionally, have included the 
cosh-like potential which serves as a good model of cold 
dark matter, known as SFCDM [3a - l37| . 

The recipe used by us to build noncommutativity into 
the FRW (flat) cosmological model was based in the ap- 
proach of reference [|[ (see also [H, HH). It can be 
summarized in the following steps: i) the Hamiltonian 
is derived from the Einstein-Hilbert action (plus a self- 
interacting scalar field action) for a FRW spacetime with 
flat spatial sections, ii) canonical quantization recipe is 
applied, i. e., the minisuperspace variables are promoted 
to operators, and the WDW equation is written in terms 
of these variables, iii) noncommutativity in minisuper- 
space variables is achieved through the replacement of the 
standard product of functions by the Moyal star product 
in the WDW equation, and, finally, iv) (semi-classical) 
modified cosmological equations are obtained by means 
of the WKB approximation applied to the (equivalent) 
modified Hamilton- Jacobi equation. 

Our results corroborate - on the grounds of a solid dy- 
namical systems basis - the results of [22[ regarding the 
dynamics at late times. Noncommutativity does actu- 
ally modify (appreciably indeed) the future asymptotics 
of the model: the late-time (future) attractor, whenever 
it exists, is associated with a solution of the cosmological 
equations that is dominated by the NC effects. Even in 
the case when the future attractor does not exist, the sta- 
bility properties of the critical points associated with the 
future asymptotics are drastically modified. This is a ro- 
bust result since it is not based on the study of a concrete 
particular solution, but on the exploration of the differ- 
ential equations flux in the phase space corresponding to 
the original cosmological model. 

However, in what regards the early-time dynamics, our 
results differ from the ones in reference |22j . Actually, 
here we have shown that the past asymptotics is also 
modified by the noncommutativity: past attractors (also, 
source critical points) are erased by the NC effects. This 
is not an unexpected result since NC effects of quantum 
nature are designed to modify the cosmological evolu- 
tion at early times, when, probably, the quantum effects 
played a major role in the gravitational dynamics. The 
discrepancy with the results of the mentioned reference 
may be a consequence of the different approaches under- 
taken. Besides, the fact that, according to the outcomes 
of [22j the early-times dynamics is not affected by the 
noncommutativity, supports our argument that the solu- 
tion studied therein was a particular, structurally unsta- 
ble solution of the modified cosmological equations. 

We want to recall that the present approach suffers 
from several "drawbacks": i) we have considered non- 
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commutativity of minisuperspace variables rather than of 
spacetime coordinates themselves, and, ii) following the 
approach of [J] we did not consider noncommutativity 
among the momenta conjugated of the minisuperspace 
variables. While the former drawback is not worrying 
since the kind of noncommutativity between the metric 
and the scalar held we have explored is expected to be 
a derived consequence of direct spacetime noncommuta- 
tivity mp, 

the latter drawback is of more concern. 
A more complete study along the lines followed here, 
where noncommutativity of the conjugated momenta is 
also considered, is the subject of ongoing research. 

Since the modified cosmological model studied in the 
present paper can be, at most, a useful toy model to seek 
for qualitative aspects of the impact of noncommutativity 
at large scales, consequently we have not discussed any 
possible observational test to check it. Anyway, given 
that inclusion of noncommutativity the way it was in- 
cluded here, affects not only the late-time cosmic dy- 
namics, but also the dynamics at very early times in a 
non-trivial way, it should be expected that the evolution 
of density perturbations in our model (in particular the 
grow of structure) is very different from the one predicted 
by other competing cosmological models, as for instance, 
the ACDM model, so that the new features can be de- 
tected by CMB measurements. The study of such an 
important issue deserves an independent publication. In 
the last instance, the present study can be considered as 
a hrst (modest) step towards a deeper understanding of 
the possible influence of primordial quantum processes in 
phenomena taking place at cosmological scales. 
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VI. APPENDIX: PHENOMENOLOGICAL TOY 
MODEL 

Here we will speculate that, instead of the scalar 
field, one has a barotropic fluid with vanishing pressure 
Pdm = (dust), that can be identified with the cold 
dark matter in the universe. To do this in a consistent 
way within the framework of our NC model, it will suffice 
to notice that for a scalar field with vanishing pressure 
p v = ip 2 jl — V = =>■ p v — (p 2 . Then it is legitimate to 
make the replacement p v — > pdm =>■ <p — > ±^/Pdm (in 
what follows, for definiteness, we shall consider only the 
positive root in the latter expression). This procedure 
will result in the following set of cosmological equations 



for our NC phenomenological model (compare with equa- 
tions (p4l) 1: 

3a 2 = p D M + Pe, 2a + 3a 2 = -p g , 

PDM + 3ap D M = -y/PDMPO, P8 = V 'PDMP9-(36) 

The first thing we want to notice is that, as before, 
the cosmic dynamics is governed by the - additional, 
non-gravitational - interaction between the CDM and 
the NC fluid. A simple inspection of the above equa- 
tions reveals that the pace of the expansion is determined 
by the correlation between the energy densities of both 
components of the cosmic mixture. Actually, it can be 
demonstrated that, in the present case, the deceleration 
parameter q = — (1 + a /a 2 ) can be expressed in the fol- 
lowing form: 

Pe — Pdm /2 

q = 1 ■ 

pe + pdm 

Hence, if the CDM energy density dilutes with the cos- 
mic expansion at higher rate than the NC fluid, as long 
as pe > pdm /2 the expansion transits from being decel- 
erated in the past to being accelerated into the future. 
The above possibility, however, depends on the way pdm 
and pg evolve with the expansion of the universe. To 
make the discussion more precise it will be mandatory 
to go onto the phase space. The asymptotic structure 
will then clearly show which solutions are generic and, 
besides, will also reveal their stability properties. 

To obtain an autonomous ODE out of the latter set 
of cosmological equations it will suffice to define a single 
phase space variable: x = ^fpg j \[3a. In terms of this 
variable the Friedmann constraint can be written in the 
following compact form: Qdm = Pdm /3a 2 = 1 — x 2 . 
Since the dimensionless CDM energy density parame- 
ter has to be necessarily a non-negative quantity, then 
| a; | < 1. However, since we are interested in cosmic 
expansion exclusively, negative values of x will not be 
considered. Another useful quantity is the deceleration 
parameter q = (1 — 3a; 2 )/2. 

The following autonomous ordinary differential equa- 
tion is obtained: 



x' = ^y xy/l - x 2 [l + v/3(l - x 2 )]. (37) 

The critical values of the x-variable are: i) x = 
- CDM-dominated (ft dm = 1), decelerated solution 
(q = 1/2), and ii) x = 1 - NC fluid-dominated (Q e = 1), 
inflationary solution (q — —1). The CDM-dominated so- 
lution is unstable. Actually, let us perturb this solution, 
i. e., x — > + e. According to (|37l) the perturbation e 
will uncontrollably grow: e(a) = e(0) exp((y / 3 + 3)a/2). 
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To explore the stability of the NC-dominated solution 
it is recommendable to make the following replacement in 
equation (|37[) : x — > y = yl — x 2 , so that the autonomous 
ODE can be written as: 2y' = -v / 3(l-y 2 )(l + V3y). We 
recall that this equation is not valid at the point y = 
(x = 1), since in the process of its derivation we divided 
by y. However, since we will be interested in perturba- 
tions around y = 0, i. e., around x = 1, but will not 
evaluate at the point y — itself, the above equation will 
be accurate enough. Lets now perform small perturba- 
tion e around y = 0. According to the latter equation 
the perturbation will decay as: e(a) cx exp(— 3a/2). In 
consequence the NC-dominated solution is stable. 



Stated in terms of the dynamical systems language: 
i) the CDM-dominated solution is the past attractor, 
while, ii) the NC-dominated solution is the future attrac- 
tor. This demonstrates that the noncommutative effects 
alone can be, indeed, a candidate to explain the late-time 
speed up of the cosmic expansion. 

The above is a nice cosmic scenario since transition 
from decelerated into accelerated expansion is generic. 
However, as with any toy model, it really does not cor- 
rectly describes the past dynamics of the universe, since 
one needs to consider, also, other cosmic components as, 
for instance, a radiation term. 



[15] 



H. S. Snyder, Phys. Rev. 71 (1947) 38. 
R. Jackiw, Nucl. Phys. Proc. Suppl. 108 (2002) 30; Phys. 
Part. Nucl. 33 (2002) S6 [arXiv:hep-th/0110057j ; L. Lan- 
dau, Zh. Eksp. Teor. Fiz. 11 (1941) 592 [English trans- 
lation: J. Phys. USSR, 5 (1941) 71]. 
C. Grojean, F. Quevedo, G. Tasinato, I. Zavala, JHEP 
08 (2001) 005 [arXiv:hep-th/0106120l . 
H. Garcia-Compean, O. Obregon, C. Ramirez, Phys. 
Rev. Lett. 88 (2002) 161301 |arXiv:hep-th/0107250l . 
H. Garcia-Compean, O. Obregon, C. Ramirez, 
M. Sabido, P hys. Rev. D 68 (2003) 045010 
|arXiv:hep-th/0210203l . 

H. Garcia-Compean, O. Obregon, C. Ramirez, 
M. Sabido, Phys . Rev. D 68 (2003) 044015 
|arXiv:hep-th/0302180"] . 

5! Estrada-Jimenez, H. Garcia-Compean, O. Obre- 
gon, C. R amirez, Phys. Rev. D 78 (2008) 124008 
|arXiv:0808.0211| . 

H. Garcia-Compean, O. Obregon, C. Ramirez, M. 
Sabido, J. Phys. Conf. Ser. 24 (2005) 203-212. 
O. Obregon, M. Sabido, E. Mena, Mod. Phys. Lett. A 
24 (2009) 1907; M. Sabido, O. Obregon, E. Mena, AIP 

Conf. Proc. 1083 (2008) 201. 

W. Guzman, M. Sabido, J. Socorro. larXiv:0812.42"5Tl 
W. Guzman, M. Sabido, J. Socorro, Phys. Rev. D 76 
(2007) 087302 [arXiv:0712.1520] . 

L.O. Pimentel, O. Obregon, Gen. Rel. Grav. 38 (2006) 
553. 

N. Seiberg, E. Wi tten, JHEP 09 (1999) 032 
[arXiv:hep-th/9908142"] . 

J. Madore, S. Schraml, P. Schupp, J. Wess, Eur. Phys. 
J. C 16 (2000) 161; B. Jurco, S. Schraml, P. Schupp, J. 
Wess, Eur. Phys. J. C 17 (2000) 521; B. Jurco, P. Schupp, 
J. Wess, Nucl. Phys. B 604 (2001) 148; J. Wess, Com- 
mun. Math. Phys. 219 (2001) 247; B. Jurco, L. Moller, S. 
Schraml, P. Schupp, J. Wess, Eur. Phys. J. C 21 (2001) 
383; X. Calmet, B. Jurco, P. Schupp, J. Wess, M. Wohlge- 
nannt, Eur. Phys. J. C 23 (2002) 363. 
A. Rezaei-Aghdam, F. Darabi, A. R. Rastkar, Phys. 
Lett. B 615 (2005) 141-145 [arXiv:gr-qc/0412089| ; 
G. D. Barbosa, Phys. Rev. D 71 (2005) 063511 



[arXiv:hep- th/040807l ]; G. Ca lcagni, Phys. Rev. D 70 
(2004) 103525 [arXiv:hep-th/0406006] ; G. D. Barbosa, 
N. Pinto-Neto, Phys. Rev. D 69 (2004) 065014; I. Hinch- 
liffe, N. Kersti ng, , Y. L. Ma, Int. J. M od. Phys. A 19 
(2004) 179-204 [arXiv:hep-ph/0205040l . 
[16] B. Vakili, N. Khos ravi, N. Khosr avi, Phys. Rev. D 
82 (2010) 103509 [arXiv: 1010. 1933] : B. Malekolkalami, 
M. Farhoudi, Cl ass. Quant. Grav. 27 (2010) 245009 
|arXiv: 1007.2499] : H. R. Sepangi, B. S hakerin, B. Vakili , 
Class. Quant. Grav. 26 (2009) 065003 |arXiv:0901.3829] : 
N. Khosravi, H. R. Sepangi, B. Vakili, Gen. Rel. Grav. 
42 (2010) 1081-1102 [arXiv: 0909. 2487] : C. Bastos, O. 
Bertolami, N. C. Dias, J. Prata, J. Phys. Conf. Ser. 174 
(2009) 012053 [arXiv:0812.3488] : O. B ertolami, C. A. D . 
Zarro, Phys. Lett. B 673 (2009) 83-89 |arXiv:0812.4607] : 
M. Maceda, A. Macias, L. O. Pimentel, Phys. Rev. 
D 78 (2008) 064041; L-H. Nguyen, R. R. Parwani , 
Gen. Rel. Grav. 41 (2009) 2543-2560 |arXiv: 0805.4294 ; 
N. Khosravi, H R. Sepangi, JCAP 0804 (2008) 011 
|arXiv:0803.T714] : L. Baro si, F. A. Brito, A . R. Queiroz, 
JCAP 0804 (2008) 005 |arXiv:0801.0810]; B. Va kili, 
Phys. Rev. D 77 (2008) 044023 [arXiv: 080 1.2438] : C. 
Bastos, O. Bertolami, N. Costa Dias, J. Nuno Prata, 
Phys. Rev. D 78 (2008) 023516 [arXiv:0712.4T22] : A. 
Bina, K. Atazadeh, S. Jalalzadeh, Int. J. Theor. Phys. 
47 (2008) 1354-1362 [arXiv:0709.3623] : M. Aguero, J. 

A. Aguilar S., C. Ortiz, M. Sabido, J. Socorro, Int. J. 
Theor. Phys. 46 (2007) 2928-2934 [arXiv:gr-qc/0703151| ; 

B. Vakili, N. Khosravi, H. R. Sepangi , Class . Quant. 
Grav. 24 (2007) 931-949 [arXiv:gr-qc/0701075 



[17] N. Khosravi, H. R. Sepangi, M. M. Sheikh- Jabbari, Phys. 

Lett. B 647 (2007) 219-224 [arXiv:hep-th/0611236| . 
[18] A. H. Cha mseddine, Phys. Lett. B 504 (2001) 33 

[arXiv:hep-th/0009153| . 
[19] J. B. Hartle, S. W. Hawking, Phys. Rev. D 28 (1983) 

2960. 

[20] J. J. Halliwell, "Introductory Lectures on Quantum Cos- 
mology" , larXiv:0909.2566l 

[21] J. Gamboa, M Loewe, J. C. Rojas, Phys. Rev. D 64 
(2001) 067901 [arXiv:hep-th /0010220| ; M. Chaichian, M. 
M. Sheikh- Jabbari, A. Tureanu, Phys. Rev. Lett. 86 



17 



(2001) 2716 [arXiv:hep-th/0010175| . 
[22] B. Vakili, P . Pedram, S. Jal alzadeh, Phys. Lett. B 687 

(2010) 119 |arXiv:1003.rT94] . 
[23] E. J. Copeland, A. R. Liddle, D. Wa nds, Phys. Rev. D 

57 (1998) 4686 [arXiv:gr-qc/9711068] . 
[24] H. Jeffreys, Proc. Lond. Math. Soc. 23 (1924) 428; L. 

Brillouin, Compt. Rend. l'Acad. des Sc. 183 (1926) 24; 

H. A. Kramers, Zeit. Phys. 39 (1926) 828; G. Wentzel, 

Zeit. Phys. 38 (1926) 518. 
[25] E. Mena, O. Obrego n, M. Sabido, Int. J. Mod. Phys. D 

18 (2009) 95 tarXiv:gr-qc/0701097l. 
[26] L. Mezincescu, |arXiv:hep-th7 0007046 
[27] M. Ryan, Hamiltonian Cosmology, (Springer, Berlin, 

1972). 

[28] D. Big gati, L. Susskind , Phys. Rev. D 62 (2000) 066004 

[arXiv:hep-th /9908056| . 
[29] M. Buric, J. M adore, G. Z oupanos, Eur. Phys. J. C 55 

(2008) 489-498 |arXiv:0709.3"i"59] . 
[30] L. Amendola, Phys. Rev. D 62 (2000) 043511 
arXiv:astro-ph/9908023] ; W 
P. Chimento, Phys 



[31] 



arXiv:astro-ph/0105479 ; 
Jakubi, D. Pavon, Phys 
[arXiv: astro-ph /0005070"] ; 
Jakubi, D. Pavon, Phys 
[arXiv:astro-ph/0303160] . 
L. P. Chimento, A. S. Jakubi, D. Pavon, W. Zimdahl 



Zimdahl, D. Pavon, 
Lett. B 521 (2001) 133 
L. P. Chimento, A. S. 
Rev. D 62 (2000) 063508 
L. P. Chimento, A. S. 
Rev. D 67 (2003) 087302 



[32] T. Gonzalez, G. Leon, I. Quiros, Class. Quant. Grav. 23 

(2006) 3165 [arXiv:astro-ph/0702227l . 
[33] B. Ratra, P. J. E. Peebles, Phys. Rev. D 37 (1988) 
3406; L. P. C himento, A. J akubi, Int. J. Mod. Phys. D 
5 (1996) 71 [arXiv:gr-qc/9506015|; A. A. Starobinsky , 
Gravit. Cosmol. 4 (1998) 88 |arXiv:astro-ph/9811360] ; 
C. Rubano, P. Scudellaro , Gen. Rel. Grav. 34 (2002) 307 
arXiv:astro-ph/0103335 ; L. A. Urena-Lopez, T. Matos, 
Phys. Rev. D 62 (2000) 081302 [arXiv:astro-ph /0003364] ; 
O. Arias, T. Gonzalez, Y. Leyva, L Quiros, Class. Quant. 
Grav. 20 (2003) 2563 [arXiv:gr-qc/0307016j . 
[34] T. Barreiro, E. J. Cop eland, N. J. Nunes, Phy s. Rev. 
D 61 (2000) 127301 [arXiv:astro-ph/9910214] ; E. J. 
Copela nd, N. J. Nunes, F. Rosat i, Phys. Rev. D 62 (2000) 
123503 [arXiv:hep-ph/0005222| . 
[35] T. Matos, L. Arturo Urena-Lopez, Class. Quant. Grav. 
17 (2000) L75 [arXiv:astro-ph /0004332] ; T. Matos, L. 
Arturo Urena-Lopez, Ph ys. Rev. D 63 (2001) 063506 
arXiv: astro-ph /0006024| 



Phys. Rev. D 67 (2003) 083513 [arXiv: astro-p h /0303 145 



[36] V. Sahni, L. M. Wang, P hys. Rev. D 62 (2000) 103517 

[arXiv:astro-ph/9910097| . 
[37] T. Matos, J. R. Luevano, I. Quiros, L. A. Urena- 
Lopez, J. A. Vaz quez, Phys. Rev. D 80 (2009) 123521 
|arXiv:0906.0396| . 
[38] M. S. Turner, Phys. Rev. D 28 (1983) 1243. 
[39] E. J. Copeland, M. Sami, and S. Tsujika wa, Int. J. Mod. 
Phys. D 15 (2006) 1753 [arXiv:hep-th/0603057] . 



